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LENSING OPTICAL DEPTHS FOR SUBSTRUCTURE AND ISOLATED DARK MATTER HALOS 
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ABSTRACT 

Multiply-imaged quasar lenses can be used to constrain the substructure mass fraction in galaxy- 
sized dark matter halos via anomalous flux ratios in lensed images. The flux ratios, however, can be 
affected by both the substructure in the lens halo and by isolated small-mass halos along the entire 
line-of-sight to the lensed source. While lensing by dark matter clumps near the lens galaxy is more 
efficient than elsewhere, the cumulative effect of all objects along the line-of-sight could be significant. 
Here we estimate the potential contribution of isolated clumps to the substructure lensing signal using 
a simple model motivated by cosmological simulations. We find that the contribution of isolated clumps 
to the total lensing optical depth ranges from a few to tens percent, depending on assumptions and the 
particular configuration of the lens. Therefore, although the contribution of isolated clumps to the lensing 
signal is not dominant, it should not be neglected in detailed analyses of substructure lensing. For the 
currently favored ACDM model, the total calculated optical depth for lensing is high, r ~ 0.2 — 20 and 
could, therefore, naturally explain the high frequency of anomalous flux ratios in observed lenses. The 
prediction, however, is highly sensitive to the spatial distribution of substructure halos in the innermost 
regions of the lens halo, which is still very uncertain. Therefore, constraints on the properties of the 
substructure population or accurate cosmological constraints, such as the mass of the warm dark matter 
particle, are difficult if not impossible to derive at this point. 

Subject headings: gravitational lensing - cosmology: theory - dark matter 



1. INTRODUCTION 

One of the generic predictions of the Cold Dark Matter 
(CDM) paradigm is a "clumpy" distribution of matter, 
with a large number of small-mass compact dark matter 
(DM) halos, both within virialized regions of larger ha- 
los (substructure) and in the field. At the same time, 
the observed number of dwarf galaxy satellites in the Lo- 
cal Group is more than an order of magnitude smaller 
than expected (Klypin et al. 1999b; Moore et al. 1999). 
This discrepancy has motivated many theoretical studies 
of alternative models designed to reduce the abundance 
of substructure (e.g., Spergel & Steinhardt 2000; Hannes- 
tad & Scherrer 2000; Hu et al. 2000) or to suppress star 
formation in small clumps via astrophysical mechanisms, 
making them dark (Bullock et al. 2000; Benson et al. 2002; 
Somerville 2002). 

If the CDM paradigm is correct, we expect ~ 2 — 10% 
of the mass of a present-day galactic halo to be tied up 
in substructure (e.g., Klypin et al. 1999b; Ghigna et al. 
2000; Kravtsov et al. 2003). The absence of apparent opti- 
cal counterparts would then indicate that the small-mass 
DM clumps are dark. In the future, it may be possible 
to detect substructure using different methods: for exam- 
ple, indirectly via the effect of clumpy matter distribu- 
tions on tidal tails of satellite galaxies (Mayer et al. 2002) 
or directly via annihilation or other observable signatures 
(e.g., Bergstromet al. 1999; Tasitsiomi & Olinto 2002, and 
references therein). At the current time, however, gravi- 
tational lensing represents the best avenue to constrain- 
ing populations of DM clumps in galactic halos (Mao & 
Schneider 1998; Metcalf & Madau 2001; Dalai & Kochanek 
2002; Metcalf & Zhao 2002). This is because substructure 
can modify the fluxes of lensed images relative to those 
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predicted by smooth lens models. The existence of such 
anomalous fluxes in many lens systems has been known 
for some time (Mao & Schneider 1998). Recently, several 
such systems were analyzed with the goal of constrain- 
ing the properties of substructure in lens halos. Dalai & 
Kochanek (2002) carried out a statistical study of seven 
radio lenses and found that the halo mass fraction in sub- 
structure, /sub, is approximately ~ 2 — 5% (see Evans & 
Witt 2003, for an alternative interpretation). This number 
is broadly consistent with the typical substructure frac- 
tions of ~ 0.02 — 0.1 predicted by numerical simulations of 
CDM models (Kravtsov et al. 2003). Note, however, that 
predicted fractions are for the total substructure popula- 
tions within the virial radius, while lensing observations 
are sensitive only to substructure in a cylinder of radius 
~ 5 — 10 kpc around the lens center. In addition, the abun- 
dance of clumps and their mass fraction can be expected to 
be depressed in the central regions due to increased tidal 
disruption and merging. 

Dalai & Kochanek (2002) and all other authors who 
have studied substructure lensing have assumed that the 
lensing is caused only by clumps in the halo of the lens 
galaxy. CDM models, however, predict that large num- 
bers of small-mass clumps also exist in the field around 
galactic halos (e.g., Sheth & Tormen 1999; Klypin et al. 
1999b). Although the lensing efficiency of an individual 
field clump is expected to be low, a typical path length 
from the source quasar may intersect the Einstein radii of 
many field clumps, resulting in a significant optical depth 
relative to that of substructure lensing. Indeed, Keeton 
(2003) recently showed that the lensing efficiency for DM 
clumps, as a function of redshift, peaks at the redshift 
of the lens but is fairly wide (Az ~ 0.3 — 0.5) and has 
tails that extend to the redshifts of the observer and of 
the source. Thus, the cumulative effect of all small-mass 
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DM clumps along the line-of-sight to the source could be 
significant. In this paper, we use a simple analytic model 
motivated by the results of cosmological simulations to es- 
timate and compare the optical depth for lensing by sub- 
structure to that of the isolated 'field' halos. We apply the 
model to two specific strong lens systems, B1422+231 and 
PG 1115+080, with anomalous flux ratios and present our 
conclusions for these systems. 

The paper is organized as follows. In section §2, we sum- 
marize the formalism of Keeton (2003), used to calculate 
the lensing optical depth of individual DM clumps. We 
describe the model for the spatial and mass distributions 
of satellite and isolated DM halos in § 3. § 4 contains the 
main details for the two observed lens systems that we an- 
alyze. Finally, we present our results and conclusions in 
§ 5 and 6. 

2. OPTICAL DEPTH 

The goal of this study is to compare the lensing effi- 
ciency of substructure located within the lens halo to that 
of clumps along the line-of-sight and in the immediate 
vicinity (but outside) the lens. We do this by comparing 
the optical depth, a measure of the lensing probability, 
of the two populations. The formalism for computing the 
lensing optical depth for small DM clumps in a strong lens 
system was recently derived by Keeton (2003) and we refer 
the reader to this paper for further details. In this section 
we summarize the main equations used in our calculations. 

Since we are concerned with spatially small perturba- 
tions to a large smooth gravitational lens, wc assume that 
lensing cross sections can be calculated by treating the 
lens as an external field. The cross-section for lensing by 
clumps is calculated assuming that their density distri- 
butions can be described by singular isothermal spheres 
(SISs). As the effect of clump lensing can vary from neg- 
ligible to very strong, we restrict ourselves to a minimum 
threshold effect on image flux (i.e., magnification or de- 
magnification). To this end, we estimate the cross sec- 
tion by taking into account all instances where the mag- 
nification or de-magnification by a clump is larger than 
a threshold value 8. We then calculate the total optical 
depth for clump lensing by integrating the cross section 
over the clump mass function in a given patch of sky and 
dividing by the area of the patch: 

t(*;«,7) = J dzD ( z ?^ J dM^a(5;K, 7 ,(3,M), 

(1) 

where k and 7 are the convergence and shear of the image 
and D{z) is the comoving distance. (5 parameterizes the 
distance ratios between the clump and the halo, so that 
in the situation where the z\ is the redshift of the object 
closer to the observer and Z2 is the redshift of the further 
object, 

where Dij = D(zi,zj) and o and s refer to the observer 
and source redshift, respectively. 

An estimate of the cross section for a general line-of- 
sight clump is given by the effective area A of the "8- 
curve" , the curve in the source plane where all images have 
magnification \i= (l + S)fio, compared to the unperturbed 



magnification 

^ = [(l-nf-^]-\ (3) 

This estimate gives the exact cross section for configu- 
rations where the clump does not change the number of 
lensed images, and a good approximation (usually a lower 
limit) for configurations where the clump creates addi- 
tional faint, unresolved microimages. Specifically, the es- 
timated cross section a can be related to A via 
det(l - I3T) 



where the matrix T is, 



r = 



det(l - r) 



A(S]K eS ,j eS ), (4) 



k + 7 
k — 7 
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The area A is a function of 8, k, and 7, where the conver- 
gence and shear are modulated by the difference in redshift 
between the lens and the clump. The effective values of 
the convergence and shear are 

(l-/3)[«-/3(« 2 - 7 2 )] 



and 



7cff = 



(1 - /3k) 2 - (/3 7 ) 2 

(1-/3)7 
(l-/3«) 2 -(/? 7 ) 2 ' 



(6) 



(7) 



The functional form of A depends on the global parity 
of the image (see Keeton 2003, eqns. 11, 15, 17, and 18 
for specific forms). For a negative parity image, the result 
is different for positive and negative perturbations. Typ- 
ically, the optical depth for the 8 < case is much larger 
than the 8 > optical depth, so clumps tend to make 
negative-parity images dimmer. 

The dependence of the cross section on the mass of the 
clump is simple and is contained within a factor of the 
Einstein radius squared, i.e., a oc A oc b 2 (M). For halos 
with a SIS density distribution, the Einstein radius is given 
by a simple expression: 
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where cr rms is the velocity dispersion of the halo, c is the 
speed of light, and Di s and D os arc the angular diameter 
distances between the lens and the source and the observer 
and the source, respectively. Due to the relative simplic- 
ity of the SIS profile and corresponding expressions, we 
assume SIS cross-sections in our analysis. Although SIS 
profiles do not describe the density distributions of cos- 
mological halos, this assumption is reasonable for our pur- 
poses, since the density distribution of total mass (baryons 
and DM) in the centers of galaxies is close to isothermal. 
Also, at this point we are primarily interested in evaluating 
the relative contributions of field clumps and substructure 
rather than evaluating the absolute optical depth. Here 
we only note that the internal density distribution of sub- 
structure halos is still uncertain and we are currently in- 
vestigating it using numerical simulations. The analysis 
presented in this paper can be easily extended for other 
profiles. 

3. POPULATIONS OF DARK MATTER HALOS 

We model the cosmological populations of dark matter 
halos by separating halos into three distinct categories: 1) 
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substructure or halos located within the virial radius of the 
lens halo; 2) halos in the vicinity of the lens: halos that 
are located outside the virial radius but in the immedi- 
ate vicinity of the lens; and 3) isolated halos or halos that 
are located far from the lens. For the latter halos, we as- 
sume a number density given by the average cosmological 
halo mass function. We do not take into account substruc- 
ture in other halos that projects onto the lens. This is a 
reasonable assumption the strong lenses analyzed in the 
literature usually have no apparent bright host projecting 
near the lens. At the same time, in cases where a bright 
galaxy is seen near the lens, it is included in the smooth 
mass model. Below we describe in detail how each halo 
category was modeled. Throughout the paper, we assume 
the currently favored flat ACDM cosmology with the fol- 
lowing parameters: h = 0.65, fiM=0.35, fiA=0.65, and 
the power spectrum normalization of er 8 =0.92 (the rms 
variance of mass distribution on the scale of 8/i _1 Mpc). 

3.1. Substructure Halos 

The substructure halo population is assumed to have a 
power-law mass function, 
dn 



dm 



(9) 



with a in the range of 1.7 — 1.9, as measured in high- 
resolution cosmological simulations (Ghigna et al. 2000). 
The spatial distribution of the substructure halos is as- 
sumed to follow that of dark matter: i.e., their spherically 
averaged number density profile is assumed to be described 
by a NFW profile (Navarro et al. 1997) outside of a core 
radius r c and constant or zero inside the core radius: 



(10) 



— oc { ( r l r ^ 1 ( 1 + r/r s ) \ r>r c 
dm \ p(r c ) or 0, r < r c . 

This model is based on the results of high-resolution cos- 
mological simulations (Colin et al. 1999; Ghigna et al. 
2000; Kravtsov et al. 2003). However, we should note that 
predictions for the spatial distribution of substructure in 
the innermost (< 0.1i? v ir) regions of halos may still suf- 
fer from "overmerging" (see discussion in § 6). As we will 
show below, our results are very sensitive to the assump- 
tions about the radial distribution of substructure halos 
and this aspect of the model is one of its largest uncer- 
tainties. 

The characteristic inner radius is given by r s = R V i r /c V i r , 
where R v i r is the virial radius of the lens halo correspond- 
ing to an overdensity of A = 340 and c„j r is the concen- 
tration parameter as given by Bullock et al. (2001), as a 
function of mass and redshift 

9 ( M \ ~°' 13 

where M* is defined as er(M*) = S c = 1.686, and zi is red- 
shift of the lens halo. The relation (9) is normalized by 
assuming that a certain fraction, / su b, of the virial mass of 
the lens is associated with substructure and setting min- 
imum and maximum limits for substructure halo masses. 
The normalization constant, K , for the relation is depen- 
dent upon the mass limits, such that 

K cx [M max 2 - Q - M min 2 ^ Q ] . (12) 

The minimum mass, M m j n , is set to be very small and has 
a negligible effect on our results for mass function slopes 



of a < 2 (it would, of course, be important for a = 2). 
For slopes that approximate the mass functions of sub- 
structure in high-resolution simulations, a 1.7—1.9, 
the normalization is sensitive only to the maximum sub- 
structure mass, M max . For our fiducial model (§ 3.4) we 
assume that M max is equal to 0.01 of the host mass. In 
cosmological simulations, M max is typically found to be 
~ 0.01 — 0.05. Satellites of larger mass are expected to 
merge with the host very quickly after accretion due to 
dynamical friction. In addition, in lens modelling the most 
massive satellites are usually incorporated as part of the 
smooth mass models. We explore the sensitivity of the 
results to M max below (see Fig. 2). 

As satellites accrete onto the lens galaxy, the matter in 
their outer portions is tidally stripped. We approximate 
the tidal radius as the radius at which the density of the 
clump equals the local density of the lens halo. The veloc- 
ity dispersions are assumed to remain unchanged. Thus, 
the mass, M, for the substructure halos in eq. 9 is the 
mass within the tidal radius. This is taken into account in 
calculating the Einstein radius and determining the mass 
limits of integration. 

Figure 1 compares the number density profile of sub- 
structure halos predicted in our model to the number den- 
sity profile measured in a high-resolution cosmological sim- 
ulation of a galaxy-sized DM halo (one of the three ha- 
los analyzed in Klypin et al. 2001). The halo was simu- 
lated in the flat ACDM cosmology (O ro = 0.3, h = 0.7, 
as = 0.9) and the profile measured at z = 0.176, typi- 
cal for lens halos. At this epoch, the mass of the halo is 
Mi8o = 1.34x 10 12 /i _1 M and the most massive substruc- 
ture halo has a mass of M max = 2.76 x lO 10 ^" 1 M Q . We 
include only halos with M > M min = 2.0 x lO 7 ^ 1 M Q , 
which corresponds to the completeness limit of the halo 
catalogs in the simulation. The mass fraction in sub- 
structure halos within the virial radius for this halo is 
/ sub - 0.087. 

The figure shows that the agreement between our model 
and the simulation result is good at large radii. At small 
radii, however, the density of substructure halos is over- 
estimated in the fiducial model; the simulation profile is 
much better approximated by a model with a larger value 
of r c . However, the simulation may actually underestimate 
the abundance of substructure in the inner regions due to 
limited mass and force resolution (see also § 6). We have 
chosen the fiducial model conservatively, possibly overesti- 
mating the effects of substructure. We estimate the effect 
of different assumptions about the number density profile 
of substructure in § 5. If it turns out that substructure 
abundance has been overestimated our conclusions about 
importance of isolated halos will be strengthened. 

3.2. Halos around the lens 

The number density profile of halos outside the virial 
radius but in the immediate vicinity of the lens are ac- 
counted for in the calculation using the cross-correlation 
function 

n m (r) = n m [l + b Ml b m £d m (r)], (13) 

where subscript m indicates the number density of objects 
of mass m, 6m, and b m are the bias for the lens halo and 
clump, respectively, and fi m is the average cosmological 
density of halos of mass m (eq. 15). We use expressions 
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equation, (1), then becomes 



Fig. 1. — The number density profile of substructure halos pre- 
dicted in our model (§ 3.1) compared to the number density profile 
measured in a high-resolution cosmological simulation of a galaxy- 
size DM halo. The halo was simulated in the flat ACDM cosmology 
(fi m = 0.3, h = 0.7, erg = 0.9) and the profile measured at redshift 
z = 0.176, typical for lens halos. At this epoch, the mass of the halo 
is Migo = 1-34 X lO 12 /*" 1 Mq and the most massive substructure 



halo has a mass of M max = 2.76 x lO 10 ^ -1 M . We include only 



halos with M n 



> M > M„ 



2.0 X 10 7 h~ l M , which corre- 



sponds to the completeness limit of the halo catalogs. The figure 
shows that the agreement between our model and simulation result 
is good at large radii. At small radii, however, the density of sub- 
structure halos is overestimated in the fiducial model; the simulation 
profile is much better approximated by a model with a larger value 
of r c . 



for the bias and the mass function from Shcth & Tormen 
(1999) (see § 3.3 for details) and the fitting formula of 
Peacock & Dodds (1996) to calculate the nonlinear dark 
matter correlation function S,dm{r)- Masses outside of the 
lens virial radius are assumed to correspond to masses of 
overdensity A = 178 as in Sheth & Tormen (1999). The 
above expression is integrated for the mass limits assumed 
in the calculation (see § 3.4). The total number density of 
dark matter halos as a function of radius n(r) is then the 
maximum of the number density from the lens substruc- 
ture model (eqs. 9 and 10) and the number density given 
by eq. 13 (see Fig. 3). 

We should note that the adopted model may overesti- 
mate the optical depth of halos within and in the vicinity 
of isolated galaxy-mass lenses. The equation 13 estimates 
the density profile around an average halo, including halos 
in very dense regions such as clusters. Thus, the aver- 
age profile n m {r) may be higher than that of a relatively 
isolated halo. Since many lenses arc found in groups and 
clusters (Keeton et al. 2000), and the two specific lenses we 
focus on here are located in galaxy groups, n m (r) should 
provide a representative profile. 

We find the optical depth by integrating the volume el- 
ement for the satellite halos - both substructure halos and 
nearby halos - over a cylinder in Euclidean geometry, cen- 
tered around the position of the image, integrating through 
the halo out to a radius of lOh^ 1 Mpc. The optical depth 



= D °-hJ 



dn 



dM—(M,x)a(M). (14) 



3.3. Isolated Halos 

For the population of isolated halos we assume a uni- 
form spatial distribution with a mass function which rea- 
sonably fits both the low mass and high mass ends of the 
mass function of halos identified in numerical simulations 
of CDM models. Specifically, we use the analytic mass 
function of Sheth & Tormen (1999) 



dn 



d lncr 



d\nM M d\nM 



where 



f {a)=A J- 
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a<5 2 



6 C 

— exp 
a 



(15) 



7^ . (16) 
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,4=0.3222, a=0.707, p=0.3, and 5 C =1.686. We use the fit- 
ting formula for the transfer function provided by Eisen- 
stein & Hu (1999) and a spherical top hat window function 
for the mass variance. As we noted above, all the halos 
are modeled as isothermal spheres, where the masses of 
the clumps are defined as in Sheth & Tormen (1999) to 
be the mass within the radius at which the overdensity is 
A=178. 

To get the number of projected isolated clumps, we inte- 
grate over the comoving volume from z = to the redshift 
of the source and over the assumed mass limits for small- 
mass clumps: 

T= r^ (z)2 ^/ dM ^ (M ' z)(7(M ' z) - 

3.4. Fiducial Model 

In our fiducial model, we set the magnification pertur- 
bation to 5 = 0.2, i.e., calculating the cross section for 
perturbations greater than 20%. The mass limits for ha- 
los outside the virial radius of the lenses are assumed to 
be 10~ 10 x M vir to 0.1 x M vir , where M vir denotes the 
virial mass of the lens. Given that halos in a lens galaxy 
which pass close to the lens center are typically stripped 
of <~ 90% of their mass by the tidal field of the host after 
1-2 orbits (e.g., Klypin et al. 1999a; Hayashi et al. 2002), 
the upper limit on substructure masses is assumed to be 
0.1 of that for the isolated clumps, i.e., 0.01 x M vir . The 
core radius of the substructure number density profile is 
assumed to be r c = 10 kpc with a constant number den- 
sity interior to it, the slope of the mass function is set to 
a = 1.8, and the fraction of mass in substructure is as- 
sumed to be / S ub = 0.1. The value of / su b is close to the 
upper end of the range found in cosmological simulations 
and maximizes the normalization of the mass function and 
the optical depth due to substructure. The sensitivity of 
the results to M max , / su b, and a is shown in Figure 2 and 
will be discussed in § 6. 

4. LENS PARAMETERS 

We estimate lensing cross-sections for two specific cases 
of quadruple- image gravitational lenses in which flux anoma- 
lies have been detected: B1422+231 and PG 1115+080. 
The source B1422+231 is a radio-loud quasar at z=3.62 
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Table 2 

Lensing optical depth due to different halo populations 
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Table 1 
Lens Parameters 



Lons/Image 
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M 


image position 
(") 


B1422+231 










A 


0.384 


0.476 


6.57 


1.014 


B 


0.471 


0.634 


-8.26 


0.961 


C 


0.364 


0.414 


4.29 


1.056 


D 


1.863 


2.025 


-0.30 


0.284 


PG 1115+080 










Ai 


0.532 


0.412 


19.96 


1.173 


A 2 


0.551 


0.504 


-19.10 


1.120 


B 


0.663 


0.644 


-3.32 


0.950 


C 


0.469 


0.286 


5.00 


1.397 



Note - k, 7, [x are the convergence, shear, and parity of the 
images, respectively, from macromodels. B1422+231 is mod- 
eled as a singular isothermal ellipsoid (SIE) plus an external 
shear (Keeton 2001). PG 1115+080 is modeled as a SIE plus 
an additional SIS representing the surrounding poor group of 
galaxies (Impey et al. 1998). 



lensed by an early-type galaxy in a poor group of galax- 
ies at 2=0.34 (Patnaik et al. 1992; Kundic et al. 1997b; 
Tonry 1998). Images A and C are bright positive-parity 
images, parity corresponding to the sign of the magnifica- 
tion, while image B is a bright negative-parity image and 
D is a faint negative-parity image (Patnaik et al. 1999). 
Keeton (2001) and Bradac et al. (2002) recently argued 
that a clump projecting in front of the image A can fit the 
lens data, concluding that mass of the perturber should be 
of order - 1O 5 ~ 6 M if it is a point mass, or - 10 6 " 7 M Q if 
the mass distribution is extended (i.e., a halo) with a SIS 
density distribution. 

PG 1115+080 is a radio-quiet quasar at z=1.72 lensed 
by an early-type galaxy in a poor group of galaxies at 
z=0.31 (Weymann et al. 1980; Kundic et al. 1997a; Tonry 



1998). Images Ai and C are positive-parity images, while 
A 2 and D are negative-parity images (Impey et al. 1998). 
Smooth lens models for PG 1115+080 arc able to fit all 
but the A1/A2 flux ratio. Table 1 summarizes the prop- 
erties of the observed images in these studies relevant to 
our study. The concentration parameter for B1422+231 
is 7.9, while that of PG 1115+080 is 7.3. The concentra- 
tion parameters are estimated using eq. 11 for the virial 
mass of the lens and its redshift. The mass, in turn, is de- 
rived from the lens model. Namely, we use the estimate of 
the Einstein radius from the lens velocity dispersion and 
invert eq. 8. The virial mass is then calculated assum- 
ing an isothermal profile as the mass within the radius 
corresponding to overdensity of A=340. The estimates of 
the virial mass and the Einstein radius of B1422+231 is 
5.1 x 10 12 M Q and 0.764" (Keeton 2001), while those of 
PG 1115+080 are 1.1 x 1O 13 M and 1.147" (Impey et al. 
1998). 

5. RESULTS 

Using the model described in the previous sections, we 
calculate and compare the lensing optical depths due to the 
substructure halos, r su b; due to all halos within r < lO/i -1 
Mpc of the lens, r<io (this optical depth includes both 
substructure and halos in the vicinity of the lens); and due 
to all 'field' halos along the line of sight from the source 
to the observer but excluding halos within 10ft" 1 Mpc of 
the lens, r>io- 

Table 2 compares r su b, t<io, and t>io for our fiducial 
model. The table shows that in the fiducial case the sub- 
structure halos provide a dominant contribution to the 
lensing optical depth for all images. The contribution of 
nearby clumps outside the virial radius of the lens (i.e., 
t<i — r su b) is negligible. The optical depth due to iso- 
lated clumps along the line of sight, r>io, is generally 
small compared to the substructure optical depth (a few 
percent) for all positive parity images. In negative par- 
ity images, the relative effect for magnification is much 
larger. Magnification, however, is a much smaller effect 
than de-magnification in these images. Thus, the contri- 
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Fig. 2. — The effect of varying a, the upper mass limit M max , and 
the core radius, r c , on the optical depth of satellites is shown here 
from left to right for B1422+231, image A, where the optical depth 
is normalized to r=l for the fiducial values. In each case, all other 
parameters of the model are kept fixed at their fiducial values. 



bution of isolated halos is relatively small. It is, however, 
not negligible and should be included in detailed analyses 
of substructure lensing. 

The results of the full numerical model can be recovered 
using simple analytic estimates. The substructure optical 
depth can be estimated as (c.f. eq. 14) 

:(Afeff) 



Tsub 



DlAx 



AM a(M eS ), (18) 



where AM w M e fj. Given that the substructure opti- 
cal depth ought to be of order unity, we can use this 
equation to estimate the effective mass. For B1422+231, 
M e ff ~ 10 8 M Q , where M c g is a tidally truncated mass. 
The corresponding virial mass of isolated halos would then 
be ~ 10 9 M Q . We can thus estimate the relative im- 
portance of the line-of-sight optical depth (c.f. cq. 17): 
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D \. Rewriting the equation explicitly in 
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For B1422+231 and the given effective mass of M e ff 
10 9 M Q , D ol ~ 10 3 Mpc and ^(M cff ) ~ IO" 9 , 

10~ 4 . The optical depth is r ~ M^ 2 , 
10 1 1 Mq the optical depth is r « 10~ 3 , roughly consis- 
tent with our numerical results. The total optical depth 
is a sum of the contributions by halos of different masses. 
This exercise shows, however, that most of the signal for 
the model we adopted is due to the high-mass end of the 
halo mass function. 

The relative contribution of isolated halos could be larger, 
if the fiducial model overestimates the abundance of sub- 
structure halos. Below we discuss the sensitivity of the 
optical depth calculations to our model parameters. The 
optical depth r su b depends on the assumed slope of the 
halo mass function, a; the mass limits, M m - m and M max ; 
the mass fraction in substructure, / SUD ; and the assumed 
shape of the radial number density profile. The lensing 
cross section a is proportional to the square of the Ein- 
stein radius, which, for SIS halos, implies 
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Fig. 3. — The number density profile for the image A of lens 
B1422+231. The two lines show the profile for substructure halos 
(eq. 10, solid line) with fiducial model parameters and the profile 
for the halos in the vicinity of the lens (eq. 13, dashed line). The 
break in the profile of the substructure halos at f» 200/i — 1 kpc (the 
virial radius) is due to change of mass from tidally truncated to 
virial (for halos outside the virial radius of the lens) in our model. 
The corresponding optical depth profile is inset. 
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Since the mass function of small-mass DM clumps is ap- 
proximated by a power-law n oc mT a , we have 

oc m 4 / 3 -, r oc [M max 7 / 3 ^ - M min 7 / 3 " 
dm L 

For a = 1.7 — 2, 7/3 — a is in the range 2/3 — 1/3, and 

the value of M m i n is unimportant. For the fiducial case, 

a = 1.8, r su b oc Mm4x 5 - The dependence on / SUD is simple: 

7*sub /"sub- 

The constraints on substructure masses for B1422+231 
suggest that the fluxes of this lens can only be affected by 
masses, M > 10 6 M Q (Keeton 2001; Bradac et al. 2002). 
Thus, masses below this value but above the minimum 
mass, should have little effect on the optical depth. The 
substructure optical depth for image A of the system for 
masses between M m ; n and 1O 6 M is t sud = 0.005, or 0.3% 
of the total optical depth. The optical depth for masses 
between M m ; n and 10 7 M Q is t sud = 0.02, or 1% of the 
total, calculated optical depth. 

The fiducial value of / SUD = 0.1 is typical for CDM ha- 
los. Although / su b varies from halo to halo, variations of 
more than a factor of ~ 2 — 3 from the fiducial model are 
expected to be rare (Kravtsov et al. 2003). Nevertheless, 
in the case when / su b and M max are at the lower end of 
the range expected for CDM halos the optical depth can 
be reduced by a factor of ~ 2. Note that these parameters 
are not independent. For example, for a < 2 lower M max 
will correspond to lower / su b- 

Figure 2 shows the effect of varying a, M max , and the 
core radius, r c , on the optical depth due to nearby and sub- 
structure halos relative to the fiducial model. The figure 
shows that varying the mass function slope a has relatively 
small effect: less than a factor of 2 for a realistic range of 
values. The variation of M max (for a = 1.8) results in only 
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a factor of ~ 2 — 3 change in optical depth. In the case 
where we model the substructure number density profile 
as described in § 3.1, our only variable parameter is r c , 
which also varies the optical depth by a factor of ~ 2 — 3. 
Thus, our results are not expected to change drastically 
due to variations of a, r c , and M max . The optical depth 
profile, however, is also affected, which we discuss below, 
and the shape of the radial number density profile is not 
well-constrained in cosmological simulations. 

Figure 3 shows the number density profile of clumps 
inside and in the immediate vicinity of the halo of lens 
B1422+231; the inset shows the cumulative optical depth 
profile for image A, which is likely to have substructure af- 
fecting its flux. The profile is normalized to unity and does 
not include the effect of isolated halos outside lOh^ 1 Mpc. 
The figure shows that most of the optical depth is con- 
tributed by regions of the highest number density. In this 
particular case, the clumps outside ~10% of the virial ra- 
dius contribute < 2% of the optical depth. Clearly, the 
substructure optical depth is very sensitive to the distri- 
bution of DM clumps in the innermost regions of the lens 
halo. In our model the inner number density of substruc- 
ture clumps is determined by the core radius r c . 

Figure 4 shows the optical depth profiles for different 
values of r c . The main panel shows the optical depth 
r(r) normalized to unity. Here again we see that halos 
at r < r c contribute most of the signal. The inset panel 
shows the corresponding unnormalizcd profiles r(r) and 
demonstrates that the total optical depth due to substruc- 
ture is quite sensitive to the inner distribution of clumps. 
The lowest curve corresponds to the case where no clumps 
are present within the central 20 kpc. In this case the op- 
tical depth is a factor of ten smaller than the optical depth 
for the fiducial model, which would increase the relative 
contribution of isolated halos to ~20-30% (c.f. Table 2). 

6. DISCUSSION AND CONCLUSIONS 

Substructure in galaxy-sized lens halos, if present, may 
alter the flux ratios of gravitationally lensed images. Stud- 
ies of multiply-imaged quasar lenses arc therefore a unique 
probe into the small-scale matter distribution in DM ha- 
los and may prove to be one of the most powerful tests 
of the CDM paradigm. Indeed, two of the most glaring 
problems for the CDM models are the density distribution 
in galaxies (predicted halos are too dense and possibly too 
cuspy) and clusters and the overabundance of small-mass 
DM clumps. Both problems are a manifestation of the 
relatively high amplitude of the CDM power spectrum on 
small scales. If DM clumps can be proven to be a unique 
explanation for the anomalous flux ratios, this would con- 
firm that the small-scale power is high and would lend 
support to the view that the problem of density distribu- 
tion in galaxies has an "astrophysical" solution. 

In order to extract useful constraints from lensing obser- 
vations, we must understand what parameters the lensing 
is sensitive to and what type of halos may contribute to the 
total lensing signal. In this paper, we presented results of 
the first study comparing the lensing optical depths due to 
small-mass halos within the lens halo (substructure), the 
halos in the immediate vicinity of the lens, and the over- 
all cosmological population of small-mass halos distributed 
throughout the entire line of sight between the source and 
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Fig. 4. — The optical depth profile for image A of lens B1422+231 
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the observer. 

The main result of our paper is that the dominant con- 
tribution to the total lensing optical depth is provided by 
the DM clumps in the densest, innermost regions of the 
lens halo. The optical depth is therefore very sensitive to 
the spatial distribution of substructure clumps within the 
lens. While the effect of masses near the lens galaxy is 
more important than elsewhere, the cumulative effect of 
all the (isolated) halos along the linc-of-sight to the source 
is not negligible. In the ACDM universe, we find that the 
contribution of isolated halos can be a sizeable fraction of 
the total lensing signal. 

The exact percentage of the isolated halo contribution 
depends sensitively on the mass function of substructure 
halos and, especially, on their spatial distribution. The 
latter is currently uncertain in cosmological simulations. 
In the highest resolution cosmological TV-body simulations 
reported to date, no clumps are typically found within the 
central 10% of the virial radius, while at larger radii the 
number density profile of substructure halos has the shape 
similar to that of the overall DM profile (Colin et al. 1999; 
Ghigna et al. 2000; Stochr 2002; Kravtsov et al. 2003). If 
such a distribution is assumed in our model, the contribu- 
tion of isolated halos to the optical depth can be as high 
as ~ 20 — 30%. The degree to which the results of cosmo- 
logical simulations are affected by resolution is currently 
unclear. It is possible that in the innermost regions the 
simulations still suffer from the perennial "overmerging" 
problem. On the other hand, the tidal disruption of NFW 
halos by the host is expected (Klypin et al. 1999a; Hayashi 
et al. 2002). In addition, clumps in real lenses would expe- 
rience enhanced tidal forces due to the baryonic material 
in the center of the lens which would make their destruc- 
tion more efficient than in iV-body simulations, and so the 
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lack of halos in the central 10% of the virial radius may 
be a real effect. 

Our results imply that flux ratio anomalies, if caused by 
DM clumps, do indeed probe the substructure population 
of lens halos and are a very promising test of CDM predic- 
tions on small-scales. For example, the alternative scenar- 
ios proposed to remedy perceived CDM problems, such 
as the simplest variants 2 of self-interacting DM (SIDM 
Spergel & Steinhardt 2000) and warm DM (WDM, e.g. 
Colin et al. 2000; Bode et al. 2001) predict a very reduced 
abundance of substructure halos and, thus, a lower lensing 
optical depth compared to that from the CDM models. 

Our estimates for the optical depth in the ACDM halos 
are in the range of ~ 0.2 — 20 and could naturally ex- 
plain the high-frequency of anomalous flux ratio images. 
In contrast, SIDM and WDM models would have opti- 
cal depths of t < 0.1, although detailed calculations are 
needed. In addition in WDM and SIDM scenarios, lensing 
by isolated halos and by substructure would have compa- 
rable optical depths because we do not expect a significant 
change in the abundance of isolated halos in these mod- 
els. To be precise, WDM models predict that abundances 
should be suppressed below a mass scale that corresponds 
to the scale of the power spectrum cut-off. The suppres- 
sion is applicable for both isolated and substructure halos. 
However, for substructure halos the abundance is further 
decreased by tidal disruption, while the abundance of the 
field halos may be increased due to fragmentation (e.g., 
Colin et al. 2000; Bode et al. 2001; Knebe et al. 2003). So 
we expect that the contribution of isolated halos would be 
even more important in the WDM models than in ACDM. 

We should also note that our results imply that precise 
cosmological constraints, such as estimates of the WDM 
particle mass or of the cross-section of the SIDM interac- 
tion, are difficult if not impossible to derive at this point. 
As shown in our analysis, isolated field halos can con- 
tribute a sizeable fraction to the lensing optical depth. 
Also, the ingredients that arc needed to estimate relative 
effects of substructure and isolated halos in various models 
are still rather uncertain. Even in the relatively well stud- 
ied ACDM cosmology, the distribution of DM clumps in 
the innermost regions of halos is not well understood. In 
particular, simulations which include the effects of baryons 
(a disk or stars in an elliptical galaxy) on the tidal disrup- 
tion of DM clumps are only starting to be employed (e.g., 
Bradac et al. 2002). Even for isolated halos the Sheth & 
Tormcn (1999) mass function that we used in our calcu- 
lations has been tested only down to M ~ 10 10 ft- _1 M Q 
(Reed et al. 2003). Given the importance of the problem, 
these issues should and will be addressed in the next gen- 
eration of cosmological simulations. Significant progress 
in determining properties of substructure halos and the 
mass function of field populations in the last several years 
makes us optimistic that uncertainties will be resolved in 
the near future. 
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2 More sophisticated SIDM models in which interaction cross-section 
varies with particle velocity may produce substructure populations 
similar to those of the CDM models (Colin ct al. 2002). 
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